Dynamical signature of a domain phase transition in a perpendicularly-magnetized 

ultrathin film 
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Domain phases in ultrathin Fe/Ni/W(110) films with perpendicular anisotropy have been studied 
using the ac magnetic susceptibility. Dynamics on time scales of minutes to hours were probed by 
quenching the system from high temperature to the stripe phase region, and varying the constant 
rate of temperature increase as the susceptibility traces were measured. The entire susceptibility 
peak is observed to relax slowly along the temperature axis, with the peak temperature increasing 
as the rate of heating is decreased. This is precisely opposite to what would happen if this slow 
relaxation was driven by changes in the domain density within the stripe phase. The data are 
instead consistent with a simple model for the removal of a significant density of pattern defects 
and curvature trapped in the quench from high temperature. A quantitative analysis confirms that 
the relaxation dynamics are consistent with the mesoscopic rearrangement of domains required to 
remove pattern defects, and that the experiment constitutes a "dynamical' observation of the phase 
transition from a high temperature, positionally disordered domain phase to the low temperature, 
ordered stripe phase. 

PACS numbers: 75.70.-i, 75.30.Kz, 75.40.Cx 



I. INTRODUCTION 

The study of domain pattern formation in two- 
dimensional systems with strong short-range attrac- 
tive interactions and weak, long-range repulsive dipole 
interactions links disparate fields such as molecu- 
lar membraneSfi crystals exhibiting high temperature 
superconductivity,— and ultrathin film magnetism.— In 
all these cases, the phase diagram is influenced by the 
strong fluctuations in two dimensions, which lead to novel 
phases and phase transitions. The domain patterns of 
ultrathin ferromagnetic films with weak perpendicular 
anisotropy have been studied intensely because this sys- 
tem is amenable to precise experimental control. Even so, 
the predicted phases and phase transitions have proven 
difficult to observe because imaging techniques average 
over the fluctuations. The purpose of the present article 
is to report observations of the slow relaxation of domain 
patterns that have been quenched from high tempera- 
ture, using the ac magnetic susceptibility. The relaxation 
is consistent with the resolution of pattern defects as the 
low temperature stripe phase is formed. 

Ultrathin films possess a surface magnetic anisotropy 
that may favour the alignment of the magnetic moments 
normal to the surface. In this demagnetization 
field can cancel much of the surface anisotropy, leaving 
a weak, temperature-dependent, residual perpendicular 
anisotropy. The resulting small domain wall energy per- 
mits magnetic domain patterns to form despite the weak- 
ness of the long-range dipole interaction. Numerous the- 
oretical analyses^i^ii and computer simulation s^'^'^°i^^ 
indicate that a "stripe" pattern is stable at low tempera- 
ture. Great progress has been made in understanding this 
phase using magnetic microscopy techniques j^'i^'i^'i^ii^ 
At higher temperature (lower anisotropy), fluctuations in 
the domain walls become important, and stripe domains 



may meander and "pinch off' to create pairs of disloca- 
tions. As the dislocations prohferate, theory predicts^i^ 
a Kosterlitz-Thouless transition to a domain phase that 
is characterized by pattern defects, domain segmenta- 
tion, and the loss of long-range positional order. At 
even higher temperature, the loss of orientational order, 
and flnally paramagnetism are predicted. As a group, 
these phases may be termed "positionally disordered" 
phases. Only a few experiments have produced magnetic 
images showing relevent defect structuresJii^ Caution 
is required in interpreting these essentially static images 
of phases driven by fluctuations. For example, images 
showing a loss of domain contrast may indicate a tran- 
sition to the paramagnetic state, or simply a dynamic 
effect that averages out the image contrast. 

A complementary method to study the domain phases 
is through their dynamics. This offers access to a wide 
range of time scales where different relaxation processes 
are active, at the expense of a pictorial representation 
of the domain pattern. There are very few studies of 
the intrinsic dynamics of these magnetic systems. Early 
work concentrated on the relaxation from the magnet- 
ically saturated statei^ Work measuring the magnetic 
susceptibility studied the transition from the stripe phase 
to a "glassy" stripe phase due to pinning of the domain 
walls by structural defectsi^i^i More recently, numerical 
simulations have suggested that the transitions between 
the various phases have distinctive dynamic signatures. 
Starting from a magnetically saturated initial state, sim- 
ulations flnd very long equilibration times for the creation 
of the stripe phase, but much shorter equilibration times 
for the positionally disordered phase. Simulations of 
quenching from a positionally disordered state indicate a 
discontinuous transition with long-lasting metastability 
before the stripe phase is establishedi^ 

The present article reports experimental studies of ul- 
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trathin Fe/Ni/W(110) films aimed at quantifying the dy- 
namics of domain pattern relaxation across a range of 
time scales by altering the constant rate of temperature 
change, R, as the magnetic susceptibility is scanned. Fol- 
lowing a controlled quenching from high temperature, the 
evolution of the susceptibility shows that the system as 
a whole slowly relaxes towards equilibrium through acti- 
vated dynamics. It is then argued that this does not rep- 
resent the relaxation of the stripe domain density (which 
occurs on a shorter time scale), but is rather due to the 
activation barriers for the mesoscopic domain wall rear- 
rangements required for the removal of dislocations and 
defects as the quenched pattern moves toward a low tem- 
perature equilibrium stripe phase. This represent a "dy- 
namical" observation of the transition from a positionally 
disordered phase to the equilibrium, stripe domain phase. 



II. THEORY 

The competition between the short range exchange 
interaction and the long range dipole interaction in a 
two dimensional magnetic system with perpendicular 
anisotropy leads to the spontaneous formation of mag- 
netic domains. A detailed study of a two dimensional sys- 
tem of domain walls, in the presence of a substrate that 
induces preferred directions of domain wall alignment, 
has outlined the expected phase diagram of the domain 
phases. ^'^'"^ At low temperature, the system forms posi- 
tionally and oricntationally ordered stripes (also termed 
the "smectic" phase) with a mean linear density UeqiT) 
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where £ is the domain wall width, Ew is the domain 
wall energy per unit length, is a constant that sets 
the scale of the dipole energy, and N is the number of 
monolayers in the thin film. The residual perpendicular 
anisotropy is reduced by spin fluctuations, so that the 
domain wall energy depends sensistively on temperature. 
This causes the domain density to increase exponentially 
with temperature. 

When a magnetic field is applied perpendicular to the 
surface, the domains with moments parallel to the field 
grow at the expense of those with moments in the an- 
tiparallel direction, creating a net moment. In the stripe 
phase, where the bending of the domain walls is un- 
favourable, the equilibrium dc susceptibility Xeq (T) is de- 
termined by the induced changes in domain wall spacing, 
and is directly related to the effective spring constant for 
stripe compression,-'?-? kefj ^ 1/Xeg, where 
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d is the film thickness. The phenomenological parameters 
In A and k represent the zeroth and first order expansion 
of the exponential dependence of the domain density on 
temperature. 



When an ac field is applied, the oscillating motion of 
the domain walls occurs through Barkhausen steps of mi- 
croscopic lengths of domain wall between thermally acti- 
vated pinning sites with time constant^! 



' pin 
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where Tpi„ is the pinning energy. Solution of a relaxation 
equation for the magnetization measured at angular fre- 
quency to gives the ac susceptibility as 



X(T)^^^^^,^Xe,(T). 



(4) 



The susceptibility falls exponentially with temperature 
on either side of a maximum. At low temperature the 
domain wall motion is stopped by pinning, and at high 
temperature by the increasing magnetic stiffness of the 
domain walls as their density increases. This characteris- 
tic shape has been observed in many studies, and permits 
a quantitative study of the pinning mechanismj^i^iiS^ 

As the temperature is increased, fluctuations in the 
domain walls are favoured, and dislocations pairs are 
formed by segmentation of the stripe domains. Accord- 
ing to the continuum model of Abanov et al^, these dis- 
locations proliferate at a Kosterlitz-Thouless transition 
to the "Ising nematic" domain phase that is character- 
ized by unbound dislocations that destroy the long range 
positional order. This (or a related positionally and rota- 
tionally disordered phase) exists for a small temperature 
range as the exponentially decreasing domain width ap- 
proaches and reaches a fundamental limit, given by the 
dipole length, or ratio of exchange and dipole energiesi^ 
Near this fundamental limit, the shape of the domain 
wall is altered2^i21 and the width becomes equivalent to 
the domain width, so that the continuum model is no 
longer strictly valid. At even higher temperature, there 
is a transition to in-plane magnetism, or paramagnetism, 
depending upon the sample thickness. 

Numerical simulations of a 2 dimensional layer of Ising 
spins are in agreement with many qualitative features of 
this description^ii^iii, particularly in observing stable 
domain patterns with the properties of the ordered stripe 
phase and of disordered phases containing defects. More 
recently, simulations have been used to study the relax- 
ation dynamics and domain phase transitions in Ising 
systems. For relaxation from a magnetically saturated 
state, to the stripe phase, Bromley et al^ find three 
distinct relaxation rates: i) for the formation of domain 
segments, ii) for the establishment of a symmetrical dis- 
tribution of up and down domains giving no net magne- 
tization, iii) for the alignment of the domains along a sin- 
gle stripe direction. Particularly the last step leads to a 
very long equilibration time. In contrast, the same study 
finds a much faster, single step relaxation to the position- 
ally and rotationally disordered phase. In simulations 
of quenching the positionally disordered phase into the 
stripe phase region, Cannas et al?- found more complex 
transitions than are predicted by the continuum model. 
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In particular, their simulations indicate a discontinuous 
transition between the stripe and nematic phases, rather 
than a Kosterlitz-Thouless transition. The discontinu- 
ous transition is accompanied by temperature hysteresis, 
and long lasting metastability of the disordered phase 
after quenching. 

Both of these numerical studies of Ising systems sug- 
gest that these distinctive dynamic responses may per- 
mit an indirect observation of different domain phases in 
experimental systems with realistic domain wall widths. 
Measurements of the magnetic susceptibility should be 
sensitive to relaxation dynamics on a number of time 
scales, but it is not certain how the magnetic susceptibil- 
ity of the stripe phase would differ from a system with 
substantial domain wall curvature and defects, such as 
the disordered domain phases. However, since the disor- 
dered phases retain short range order, and because the 
domain density is expected to increase with temperature 
irrespective of the phase, the ac susceptibility should con- 
tinue to be described, at least qualitatively, by eq.(|3]). On 
the other hand, when domain wall fluctuations, curvature 
and dislocations become important and contribute to the 
elastic energy of the system, one should expect quanti- 
tative changes in the effective spring constant kejf, and 
therefore in the susceptibility. 



III. EXPERIMENTAL RESULTS 

Films of 1.5 ML Fe/2.0 ML Ni/W(110) were grown in 
ultrahigh vacuum and studied in situ. The growth proce- 
dures were taken from previous studies of their structural 
and magnetic propertiesi^ The films have perpendicular 
anisotropy at low temperature. For films thicker than 2.2 
ML, there is a re-orientation transition from perpendicu- 
lar to in-plane moments as the temperature is increased, 
followed by a Curie transition to paramagnetismi^^ The 
film thickness, cleanliness and structure were confirmed 
using Auger electron spectroscopy and low energy elec- 
tron diffraction, ac magnetic susceptibility measure- 
ments at a frequency of 210 Hz using a field amplitude 
of 2.0 Oe were made using the surface magneto-optical 
Kerr effect, and lock-in amplification.^ - The sample tem- 
perature was measured using a W-Re5%/W-RelO% ther- 
mocouple, and was controlled by radiant heating from a 
filament, and conductive cooling through a copper braid 
attached to a liquid nitrogen reservoir. When heating, 
the average rate of change of the sample temperature, R, 
could be controlled from 0.05 K/s to 1.00 K/s, with fluc- 
tuations in the rate less than 0.05 K/s. The maximum 
rate of controlled cooling that could be maintained over 
the entire relevant temperature range was -0.10 K/s. 

A number of steps were taken to ensure reproducible 
measurements. First, the films were annealed to 400K 
after growth, and subsequent measurements did not ex- 
ceed 360 K. In each case where the susceptibility was 
measured by heating from low temperature, the sample 
was first cooled from 360 K at a rate of R = -0.10 K/s. 
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FIG. 1: Magnetic susceptibility of a 1.5 ML Fe/2.0 ML 
Ni/W(110) film measured at different constant rates of tem- 
perature change, R. a) Measurements for heating the film, 
taken after cooling from 360 K at -0.10 K/s. b) The tem- 
perature where the susceptibility peaks, as a function of the 
heating rate. The fitted line is discussed in section IV. 



Since the initial cooling of the sample and sample support 
was much slower, the first measurement upon heating was 
discarded. When a series of curves were measured from 
the same film at different heating rates, the order of the 
rates was randomized so that effects of aging could not 
masquarade as effects due to changes in the heating rates. 
Finally, the effect of changes in the time constant used 
for lock-in amplification in the experiments was studied. 
In this way the time constant was chosen as 2 s. This is 
as large as possible to reduce noise, yet smaller than the 
value that affected the shape of the experimental curve. 

The real part of the ac magnetic susceptibility mea- 
sured at a series of heating rates on a single film is shown 
in fig. la. The shape of the curve changes only in subtle 
ways as the heating rate is changed - the primary effect of 
changing the heating rate is a shift of the entire suscepti- 
bility curve in temperature. In the range 0.70 K/s > R > 
0.30 K/s, there is very little shift, but for slower heating 
rates the peak shifts progressively to higher temperature, 
and the susceptibility curve broadens somewhat. The lo- 
cation of the peak as a function of the heating rate is 
summarized in fig. lb. 
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FIG. 2: a) Representative fits of the susceptibility curves in 
fig. la to eq.Q. b) The parameter k in the fit. The line is a 
guide to the eye. c) The parameter Tpin in the fit. The line 
is the average value. 



The changes in the susceptibility curve upon heating 
have been characterized by fitting the curves to eq.([2]|4]). 
The fits are summarized in fig. 2. In part a), the solid 
lines show the fits to the data for the extreme values of 
R from fig. 1. These fits are representative, with the 
curve following the data closely except at the low tem- 
perature tail. The low temperature tail deviates from 
the data because the fit assumes a single average pinning 
energy, Tpin, rather than a more correct distribution of 
pinning energies. There are relatively few pinning sites 
with a pinning energy significantly lower than the aver- 
age, but these are the only sites where the domain walls 
remain free at low temperature, and thus account for all 
the signal in the tail of the curvei^i The greater size of 
the low frequency noise in the data for R = 0.70 K/s is 
consistent with the reduced measurement time compared 
to the data with R = 0.03 K/s. Fig. 2b and 2c present 
the fitted parameters k and the average defect pinning 
energy Tpm, respectively. It is clear that the shift in the 
peak temperature is associated with a change in k. The 
pinning energy remains almost constant, as is expected, 
since the pinning energy is a characteristic of the film 
and substrate microstructure and does not change with 
the heating rate. 



Since k is an effective parameter, it is difficult to be 
precise about what causes it to change with heating rate 
so that the peak shifts in temperature. However, a num- 
ber of possible artifacts can be ruled out. Fundamentally, 
changes in the equilibrium domain patterns (and thus 
the susceptibility) are driven by changes in the perpen- 
dicular anisotropy as a function of temperature. How- 
ever, the perpendicular anistropy is determined by ma- 
terials properties (magneto-crystalline anisotropy, shape 
anisotropy) that are not dependent on the rate of tem- 
perature change. The peak shifts are therefore not due 
to a true shift in the temperature dependence of the per- 
pendicular anisotropy. It is, however, conceivable that 
the peak shifts are due to a time lag between the mea- 
sured and true temperature, or due to the response time 
of the temperature measurement system itself. A sim- 
ple calculation of the thermal diffusion time across the 
tungsten substrate gives an upper limit of a few tenths 
of K for a shift due to the fact the temperature mea- 
surement point and experimentally probed area of the 
film are physically separated. Limits on the size of the 
thermal lag of the overall temperature measurement sys- 
tem were established by heating at a constant rate up 
to a specific temperature, and then monitoring the mea- 
sured temperature as a function of time at a constant 
temperature. In this case, a conservative upper limit in 
the temperature shift of 1 or 2 K was established. These 
effects are much smaller than those shown in fig. lb. 

It can also be shown that the changes in k with heat- 
ing rate are not due to the dynamical processes involved 
in changes to the domain density with temperature. At 
low temperature, the domain density will be small. As 
the temperature is increased, the domain density will in- 
crease through the creation and growth of new stripes. 
Since this process is thermally activated, the domain den- 
sity will always lag the equilibrium domain density. As 
the heating rate is increased, this effect will be more pro- 
nounced, and the amount of lag between the actual and 
equilibrium density at a given temperature will also in- 
crease, and this could affect the fitted value of k and 
shift the susceptibility curve. However, the crucial point 
is that, during heating, the lagging domain density is 
smaller than the equilibrium value, so the measured sus- 
ceptibility is increased, causing the measured susceptibil- 
ity curve to be progressively shifted to higher tempera- 
ture for higher heating rates. This is precisely the op- 
posite of what is seen experimentally. Although changes 
in domain density must occur, apparently this dynam- 
ical process is active on a faster time scale than those 
probed by the present experiments on these samples. On 
the time scale of the present measurements, the system 
is able to continually maintain a domain density close to 
the equilibrium value, except perhaps at the lowest tem- 
peratures. This conclusion is further corroborated by the 
very good fits of the susceptibility curves to eq. ([4]) reflect- 
ing the underlying processes of domain wall pinning at 
low temperature and a rapid increase in domain density 
at high temperature. The shift in the peak temperature 
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FIG. 3: The magnetic susceptibility measured while the sam- 
ple was cooled at a few different rates. 

in fig. 1 must be predominantly due the relaxation of 
some quantity other than domain density. 

A series of further experiments were performed to more 
completely characterize the system. The susceptibility 
measured for a different sample during cooling is pre- 
sented in fig. 3. Even though the range of R available 
for cooling is limited, it is clear that the difference in 
the curves for R = -0.10 K/s and -0.05 K/s is at most 
small and likely negligible. There is a distinct asymmetry 
in the behaviour for heating and cooling at these rates. 
The peak position and shape for the cooling curves is 
most similar to those of the heating curves at R > 0.40 
K/s. 

Longer time scales have been probed by stopping the 
heating or cooling at a predetermined temperature, and 
monitoring the susceptibility as a function of time. These 
results are presented in fig. 4a and 4b, for heating and 
cooling respectively, with \R\ = 0.10 K/s. After heating 
to 252 K, the susceptibility relaxes to a lower value by 
a simple exponential decay with Tr — 618 ± 3s, whereas 
after heating to 282 K, the susceptibility relaxes more 
quickly (r^ = 297 ± 3s) to a higher value. These results 
are in agreement with those presented in fig. la; the 
relaxation is always toward the curve measured with a 
smaller heating rate, which is closer to equilibrium. The 
curves in fig. 4b confirm the asymmetry between heat- 



FIG. 4: a) The magnetic susceptibility as a function of time, 
measured while heating at 0.10 K/s. The heating was stopped 
at the specified temperature, and then monitored as a func- 
tion of time at constant temperature. The solid line is an 
exponential fit to the relaxation, b) As in a), except that 
the susceptibility is measured while cooling at -0.10 K/s. c) 
A different sample is heated at 0.30 K/s (dashed line), then 
cooled at -0.10 K/s and measured during heating again (solid 
line). This time the heating is stopped at 282 K, and the 
susceptibility is monitored as a function of time at constant 
temperature (inset). The heating is then recommenced, pro- 
ducing the remainder of the curve. 



ing and cooUng. When coohng is stopped at these same 
temperatures, the long-term relaxation is in the same 
direction as for heating (after a transient), but the relax- 
ation time constant is much larger. Rough estimates are 
Tr w 5, 000 s at 250 K and « 1, 200 s at 281 K. 

A final experiment on a third film is presented in fig. 
4c. In this case, the susceptibility was measured for a 
heating rate of i? = 0.30 K/s (dashed curve), the sam- 
ple was cooled once more and the measurement was re- 
peated (solid curve). This time the heating was stopped 
at T=282 K, and the relaxation of the susceptibility was 
monitored as a function of time, as indicated in the in- 
set. After the relaxation was complete, heating at R = 
0.30 K/s was resumed, producing the remainder of the 
solid curve. This result confirms that the shift in the 
curves is not caused by the different heating rates per se, 
but rather that the system is relaxing from a state that 
produces a susceptibility peak at lower temperature to 
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one with a peak at a higher temperature, and that the 
amount of relaxation depends upon the total time that 
has elapsed. 



IV. ANALYSIS AND DISCUSSION 

These results provide evidence that some quantity in 
the quenched samples is relaxing as the equilibrium stripe 
phase is established, and that this relaxation occurs on a 
longer time scale than that required for the relaxation of 
the domain density. This is qualitatively consistent with 
the numerical results of Cannas et al^ who suggest the 
high temperature, positionally disordered smectic phase 
will be strongly metastable and persist into the stripe 
phase. It is also qualitatively consistent with Bromley et 
al^, who find that the resolution of an initial state with 
many pattern defects and dislocations into an ordered 
stripe pattern proceeds with a very long time constant. 
We therefore advance the hypothesis that quenching the 
system traps a configuration with significant domain cur- 
vature and density of domain dislocations at low tem- 
perature, and that the system relaxes very slowly to an 
ordered stripe pattern. Upon heating the sample, the ex- 
periments indicate a more efficient resolution of the pat- 
tern defects during domain creation than was the case 
during domain annihiliation (cooling). This more effi- 
cient relaxation is revealed by a shift of the susceptibility 
peak to higher temperature, and a change in the phe- 
nomenological parameter k, that is related to domain 
formation. It allows the system to approach the equilib- 
rium state at high temperature, so that the system shows 
reproducible and systematic hysteresis upon temperature 
cycling with different heating rates. 

To test this interpretation, it is first assumed that the 
relaxation rate during cooling is so much slower than dur- 
ing heating that the amount of relaxation that occurs 
during cooling can be neglected. The relaxation during 
heating is modeled as an activated process with time con- 
stant 



Tr = Tor e:iq>{Tr/T), 



(5) 



where Tr is an activation barrier that must be overcome 
in relaxing to the equilibrium stripe phase. The total 
number of time constants that have passed while heating 
from the initial temperature Ti up to the peak tempera- 
ture Tpeak can be defined as t^f / : 



teffiR) 



dT 



RTr{T) 



(6) 



Using the peak of the susceptibility curve as a marker to 
follow the relaxation suggests 



Tpeak{R) = To - Aexp(-te//(i?)). 



(7) 
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To is the peak temperature when relaxation is com- 
plete. A plot of ln(To — Tpeak) vs. teff has a slope of 
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FIG. 5: Fit of the relaxation of the temperature of the suscep- 
tibility peak to eq.((7)|. a) Peak temperatures plotted agains 
te// as defined in eq.®, assuming tor = 0.7 s. Tr = 1637 
K is the fitted value of the nucleation energy for the phase 
change. To = 270.5 K and A = 25.0 K. b) The residuals of 
the logarithmic fit in part a), as a function of logjQ(ror). 



-1, and intercept of In A. There are three adjustable pa- 
rameters, TQr, Tr, and Tq, but two are linked by the re- 
quirement that the slope of the plot is -1, leaving two in- 
dependent parameters. The least squares fit of the peak 
positions to eq.([7]) is shown in fig. 5a, and as the solid 
line in fig. lb. The value chosen for tq,- is constrained by 
the least squares residual in logarithmic space, as shown 
in fig. 5b, which gives a lower limit of ror ~ lO'^, but does 
not establish an upper limit. However, the two experi- 
mental relaxation times in fig. 4a allow an independent 
estimate of the two quantities that determine the time 
constant in eq.® as tq^ = 0.7 s, and Tr = 1735 K. This 
value of Tor is entirely consistent with the lower limit es- 
tablished by the fitting, and was used for the plots in fig. 
5a and lb. The fitted value of Tr = 1637^^^, correspond- 
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ing to Tor = 0.7 s, is in very reasonable agreement with 
the experimental estimate of 1735 K. It is clear that this 
simple analysis gives a very good account of the data. 

The excellent fit to the simple relaxation model for the 
peak as a whole has a number of implications. First, 
it confirms that "quenching" is an appropriate descrip- 
tion of the cooling of the system, since essentially no 
relaxation occurs during cooling. The measurement of 
x{T) for R = 0.70 K/s in fig. 2a is, in some important 
way, unrelaxed from the initial high temperature state 
before quenching. With teff — 0.07 time constants, ap- 
proximately only 7% of the relaxation has occured. By 
contrast, the measurement with R = 0.03 K/s has t^ff 
= 3.52 time constants, so that 97% of the relaxation is 
complete. The data at intermediate values of R in fig. la 
are not a linear combination of these two endpoints - if 
they were, x(T) would exhibit a double-peaked structure. 
It therefore seems that the relaxation is not related to a 
localized property such as, for instance, phase separation 
or nucleation and growth of a low-temperature phase. 
Rather, a distributed property, which is ultimately re- 
lated to K, is relaxing. The analysis is therefore consis- 
tent with the relaxation of a non-equilibrium density of 
domain dislocations, and/or domain wall curvature that 
are characteristic of high temperature state from which 
the quench occurred. 

The fundamental time scale of the relaxation, ror, is 
of order 1 s, which is very long compared to the funda- 
mental time scale for localized Barkhausen steps in the 
domain wall position, Topm ~ 10~^ s, in eq.([3]). This im- 
plies that the relaxation requires changes of the domain 
walls on a mesoscopic scale involving the co-ordinated ac- 
tion of many independent microscopic Barkhausen steps, 
and is not directly related to the shape profile of the do- 
main wallj^i^ The activation energy, Tr w 1600 K, for 
this co-ordinated motion of many Barkhausen steps, is 
about 1/4 of the activation energy Tpm ~ 6200 K for 
the pinning of just one local Barkhausen step. Taken to- 
gether, these indicate that the relaxation involves larger 
scale domain rearrangments driven by an interaction that 
is intrinsically weak but becomes substantial when inte- 
grated over a large area. This again is consistent with the 
relaxation of domain curvature and dislocations, driven 
by dipole energies. 

In addition, there is clearly a relation between the 
phenomenological parameter k and the relaxation. The 
larger values of n for susceptibility measurements with 
larger values of R indicate that after quenching, but be- 
fore relaxation, the domain pattern is magnetically stiffer 
- that is, the effective spring constant for the compression 
of the pattern by a magnetic field is significantly larger. 
In the equilibrium stripe phase at low temperature, the 
magnetic stiffness depends only upon the domain density. 
Curvature of the domain walls and dislocations are essen- 
tially absent in equilibrium because they would greatly 
increase the elastic energy. At high temperature and in 
the positionally disordered phases, on the other hand, 
the enhanced fluctuations of the domain walls can accom- 



modate curvature and an equilibrium density of domain 
dislocations without a significant cost in elastic energy. 
If quenching from high temperature introduces a rela- 
tively large amount of domain curvature and density of 
dislocations to the low temperature domain pattern, it is 
reasonable that this would be reflected in a larger elastic 
constant and value of k, which would decrease, as is ob- 
served in fig. 2b, as the relaxation to equilibrium removes 
the defects. 

Finally, there is the question of the asymmetry of 
the relaxation times for heating and cooling the system. 
Since domain annihilation upon cooling proceeds by first 
breaking up existing domains, it creates additional pat- 
tern defects and works against the removal of existing de- 
fects, making the relaxation very slow. Domain creation 
during heating, on the other hand, provides an oppor- 
tunity for defects to be removed by growth and merging 
with other defects, or at the very least dilutes them in a 
greater stripe density. This abets the relaxation of the 
system toward equilibrium. 



V. CONCLUSIONS 

Significant hysteresis and a heating/cooling asymme- 
try has been observed in the magnetic susceptibility of 
perpendicularly magnetized 1.5 ML Fe/2 ML Ni/W(110) 
films. The entire susceptibility curve relaxes slowly to 
higher temperature as the rate of heating during the mea- 
surement is decreased. This is incompatible with relax- 
ation of the domain density, which occurs on a faster 
time scale and causes a temperature shift of opposite 
sign. However there is a great deal of self-consistent ev- 
idence that it is compatible with the relaxation of dislo- 
cations, defects and curvature that have been quenched 
from high temperature into an unstable domain config- 
uration at low temperature. This evidence includes; an 
excellent quantitative fit of the relaxation to an activated 
relaxation with a fundamental time and energy scales ap- 
propriate for a large scale domain rearrangement by the 
weak dipole interaction; the correlation of the relaxation 
with a reduction in the magnetic "stiffness" of the film; 
and a qualitative link between the annihilation/creation 
of domains and the creation/annihilation of defects to 
understand the asymmetry in heating and cooling . This 
suggests that in quenching from 360 K to 220 K, the 
system passes through a transition from a positionally 
disordered domain phase to the semectic stripe phase. 

These observations give qualitative support to recent 
numerical studies that indicate that the removal of do- 
main dislocations after moving across a phase transition 
to the ordered stripe phase is very slow. However, it 
is not possible to discriminate between the detailed pic- 
tures suggested in those studies; whether the transition 
is of the Kosterlitz-Thouless type or is discontinuous^-, 
whether it originates in the Ising nematic or tetragonal 
domain phasea^^, and whether the quenched system is 
metastable or merely dynamically hindered. Regardless 
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of these details, the results have important implications 
for the interpretation of studies using imaging techniques 
to study these domain phases. They suggest that great 
care must be taken to obtain images of phases that rep- 
resent equilibrium conditions. Most imaging studies do 
not give details of the thermal history of the samples, 
although some mention very long times for the evolu- 
tion of some patterns. It is often not clear that images 
that seem to show paramagnetic or tetragonal phases are 
close to equilibrium and can serve as a basis for creating 
a phase diagram. It may be that the imaging of patterns 
with many defects is possible only because of the very 



slow dynamics observed in the present study. 

The present dynamical studies highlight two different 
processes that occur on very different time scales in per- 
pdendicularly magnetized domain patterns - the motion 
of domain walls in response to an applied field and the re- 
moval of pattern defects in response to long range dipole 
interactions. We expect that there will also be a dynami- 
cal signature of the annihilation and creation of domains 
to adjust the domain density on an intermediate time 
scale, and are currently conducting experiments to ex- 
amine and characterize this process. 
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